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ABSTRACT 

The nuclear symmetry energy figures crucially in the structure of asymmetric nuclei and, 
more importantly, in the equation of state (EoS) of compact stars. At present it is almost 
totally unknown, both experimentally and theoretically, in the density regime appropriate 
for the interior of neutron stars. Basing on a strong-coupled structure of dense baryonic 
matter encoded in the skyrmion crystal approach with a topology change and resorting to 
the notion of generalized HLS (hidden local symmetry) in hadronic interactions, we address 
a variety of hitherto unexplored issues of nuclear interactions associated with the symmetry 
energy, i.e., kaon condensation and hyperons, possible topology change in dense matter, 
nuclear tensor forces, conformal symmetry and chiral symmetry etc in the EoS of dense 
compact-star matter. One of the surprising results coming from the hidden local symmetry 
structure that is distinct from what is given by standard phenomenological approaches is 
that at high density, baryonic matter is driven by RG flow to the "dilaton-limit fixed point 
(DLFP)" constrained by "mended symmetries." We further propose how to formulate kaon 
condensation and hyperons in compact-star matter in a framework anchored on a single 
effective Lagrangian by treating hyperons as the Callan-Klebanov kaon-skyrmion bound 
states simulated on crystal lattice. This formulation suggests that hyperons can figure in 
the stellar matter - if at all - when or after kaons condense, in contrast to the standard 
phenomenological approaches where the hyperons appear as the first strangeness degree of 
freedom in matter thereby suppressing or delaying kaon condensation. In our simplified 
description of the stellar structure in terms of symmetry energies which is compatible with 
that of the 1.97 solar mass star, kaon condensation plays a role of "doorway state" to 
strange-quark matter. 
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1 Introduction 



The question as to whether strong interactions have anything to do with black hole 
formation or rather with the other side of the coin, the maximum mass of compact stars 
stable against collapse to black holes, can be addressed in terms of the nuclear symmetry 
energy egym- How nuclear forces mediating the strong interactions and controlling the phase 
structure of multi-baryonic systems enter into the nature of egyrm presumably controlling 
the fate of compacts stars, is one of the principal themes of our current theoretical research 
into baryonic matter at high density. This note is a sequel to, and a substantial updating 
of, the previous reports [1, 2]. 

On the experimental side, there are broadly two avenues, namely, terrestrial labo- 
ratories and space observatories. The heavy-ion collision is one of the terrestrial tools to 
gain information on how the strong interactions between hadrons (nuclear interaction) are 
modified with varying densities. In project are such laboratories as KoRIA (Korean Rare 
Isotope Accelerator, officially called "RAON") and other RIB machines, FAIR (Facility for 
Antiproton and Ion Reactions) at GSI/Darmstadt, FRIB(Facility for Rare Isotope Beams) 
at MSU/Michigan, and NICA (Nuclotron-based Ion Collider Facility) at JINR/Dubna. 
They will probe the ranges of temperature up to 50 — 60 MeV and of density up to (3 — 4)no 
(where no is the normal nuclear matter density no ~ 0.16 fm~^) in the nucleon number 
density [3]^ commensurate with the temperature and density encountered in binary mergers 
of neutron stars and black holes [4] within the transient time of 10 — 20 fm/see [5]. The n-p 
asymmetry, 6{= = 1 - 2x,x = P/{N + P)), ranges from for to 0.198 for ^^^Au 
or 0.227 for ^^^U collisions [6]. 

In heavy-ion collisions the transient time for the dense matter phase is too short to 
activate the weak interactions. Hence the initial n-p symmetry factor 6 does not change, 
so any process involving net strangeness number production is suppressed. This differ- 
ence between the hadronic matter in heavy ion collision and the hadronic matter in weak 
equilibrium of stellar matter renders the role of the symmetry energy different from each 
other. 

For the stellar mass less than a half of solar mass, the relevant density is less than uq. 
The range of density probed by the future facilities like FAIR and NICA is quite relevant 
for the neutron star with a mass of order of solar mass, where the central density is of order 
of a few times no [7]. The newly observed neutron star of 1.97 sola mass is consistent with 
the higher central density with the upper limit ~ lOno [8] [9]. The symmetry energy and 
the effective masses of hadrons in medium (including koan) are the quantities depending on 
how strong interactions are implemented in the dense hadronic matter. 

In this paper, we make a highly simplified discussion of how flavor symmetry can 
figure in the equation of state (EoS) that enters into the structure of compact stars, i.e., 
their maximum mass stable against gravitational collapse vs. their radii. Our central theme 

^In this paper, unless otherwise noted, baryon number density will be denoted n. 
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is that baryonic matter at densities exceeding that of normal nuclear matter at present is 
more or less unknown, and that one promising approach is to put skyrniions (in 4D) or 
instantons (in 5D) on crystal lattice [10]. We anchor our reasoning on the hidden local 
symmetric nature of baryonic interactions that has been discussed since 1980's but recently 
rediscovered in gravity-gauge dual approach to the strong interactions where an infinite 
tower of hidden gauge fields figure. 

We admit that given that the work which is technically and numerically involved is 
in an early stage, most - if not all - of what we find in this paper are far from rigorous, 
drawing ideas from qualitative or at best semi-quantitative results so far obtained in the 
WCU-Hanyang project. Even so, some are novel and promising. 

The principal results we have obtained so far are as follows. 

Our argument is largely based on the observation that when put on crystal lattice, a 
matter of multi-skyrmions (or multi-instantons in 5D) representing dense baryonic matter 
makes a topological phase change to a matter of half-skyrmions (or "dyonic salts") at a 
density ni/2 lying (not far) above the nuclear matter density uq. While this statement 
is firm only in the limit, the topology change is likely to be robust and its physical 
implication could be independent of the crystal artifacts. 

A prediction [12, 13] that follows from the hidden local symmetric nature of the strong 
interactions that is distinct from other approaches is that the compressed baryonic matter 
is driven via "symmetry mending" [14] to an infrared fixed point called "dilaton-limit fixed 
point" which lies at a density near but below chiral restoration. The consequence is that 
since the isovector channel involving p-NN coupling is affected, the standard RMF approach 
where the p-NN coupling controls the symmetry energy becomes a suspect. 

One of the striking effects of the topology change is the modification of what was 
formerly known as "BR (Brown-Rho) scaling " to a new scaling that will be referred, for 
reasons that will be made clear, to as "BLPR (Brown-Lee-Par k-Rho) scaling." This has 
a drastic consequence on the nuclear symmetry energy. Due to this topology change, the 
structure of nuclear tensor forces gets significantly altered, which then determines the struc- 
ture of the nuclear symmetry energy which is dominated by the tensor forces. The topology 
change also influences how hadrons behave on top of the dense background, such as the 
onset of kaon condensation and the appearance of hyperons at high density, a controver- 
sial issue for the EoS of compact stars. We suggest that kaon condensation and hypcron 
presence in compact-star matter can be treated on the same footing in a unified way. Our 
approach predicts that hyperons appear only after kaons condense in disagreement with the 
usual phenomenological mean field approaches where hyperons appear first as strangeness 
degree of freedom and then suppress or delay kaon condensation. 

We propose that kaon condensation can be taken as a doorway to strange quark 
matter, a notion which seems to be consistent with the structure of the 1.97 solar mass 
star. [11] 
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2 Symmetry Energy 



The nuclear symmetry energy denoted egym figures importantly in nuclear physics 
and in compact-star physics. Although it is more or less controlled by experiments up to 
near nuclear matter density no - and is probed slightly above no in heavy-ion collisions, 
it is almost completely unknown at high density relevant to the interior of compact stars 
going up to, e.g., ~ lOno- There are a large number of theoretical predictions for egym that 
range widely above no but given the total absence of model-independent tools to gauge 
the reliability and experimental information, none of them can be taken seriously. In this 
section, we discuss the various aspects of the symmetry energy, some rigorous and some 
not, that have not been addressed in the literature. 

2.1 Isospin symmetry 

To start with, we review the elementary - albeit text-book - notion of the symmetry 
energy so as to define the language that wc shall use for our analysis. Although there is 
nothing new in this discussion, it will sharpen our arguments to be developed in a way to 
enable us to address the issues raised in Introduction. 

The (strong) interaction between nucleons is charge symmetric and has been formu- 
lated in an SU{2) symmetric way, i.e., isospin symmetry. Proton and neutron are assigned 
to be members of a doublet and all hadrons are classified into SU (2) multiplets. 

Nuclei are the systems where many nucleons are bound together. Involving many- 
body correlations, it is difficult to describe them directly from the basic nucleon-nucleon 
interaction, not to mention from QCD. Even when we are provided with an effective La- 
grangian or Hamiltonian for nucleon-nucleon interactions, it is not a straightforward task 
to calculate physical observables of many-body systems from it. 

Throughout this paper, we will ignore Coulomb interactions. Then the Lagrangian is 
SU (2)-symmetric, which is blind to flavor in the sense that it cannot distinguish between 
proton and neutron. The interaction Hamiltonian commutes with SU{2) generators, /, 

[Ii,Hint]=^. (1) 

[Ii,lj]=ieijklk- (2) 

Thus the proton (I3 = 1/2) numbers and neutron (I3 = —1/2) numbers are conserved, 
so that we can classify the eigenstate by the definite number of protons and neutrons 

\Np,Nn). (3) 

The energy of the eigenstate of the Hamiltonian does not depend on 73 but on P. The 
eigenstate, Eq.(3), can be decomposed into the irreducible representations (multiplets) of 
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SU{2), Clebsch-Gordan series, as 

\Np,Nn) = ^iCi\I :Np,Nn) (4) 

where |/p = /(/ + 1). The energy of the state is given by 

E{Np,Nn) = {Np,Nn\H\Np,Nn) = ^i\Ci\^Ei (5) 

where Ej is a reduced matrix element of the Hamiltonian, H = Hq + Hint, 

Ei = {I\\H\\I) (6) 

which is independent of /s, and depends on the details of the strong interactions for each 
/-channel. 

Although it may appear that the energy is independent of the compositions of protons 
and neutrons for a given total number of nucleons 

N = Nj, + Nn, (7) 

different compositions have different decompositions into multiplctcs, i.e., different sets of 
C/. Therefore different compositions of protons and neutrons lead to different energies. 
Here are a few examples: 

1. The states with Np and Nn exchanged have the same set of C/, so the energy is 
degenerate, 

E{Np,Nn) = E{Nn,Np). (8) 

2. The state with Np = N and Nn = has the highest I3 but the state with Np = N/2 
and Nn = N/2 cannot have that state. Moreover, the state with Np = N/2 and 
Nn = N/2 can have an isosinglet component while the state with Np = N and 
Nn = cannot. Therefore they should have different decompositions of multiplets or 
equivalently different sets of Cj. Hence, in general, 

E{Np = N,Nn = 0)^ E{Np = 1/2, 7V„ = 1/2). (9) 

This is of course a well-known relation that explains why the nuclear symmetry energy 
appears in asymmetric nuclear matter. It is not a result of isospin symmetry breaking 
of the strong interactions. The strong interaction is isospin-symmetric, but we are 
considering the states with different Clebsch-Gordan decompositions. 

In short, the density dependence of nuclear symmetry energy can be understood in 
terms of the density dependence of C/(n) and Ei{n). The normalization, S/IC/p = 1, im- 
plies that at higher density, the probability of having a specific / becomes very small. How- 
ever, since the total energy of the system depends on the product of Ci{n) and {I\\Hint\\I) ■, 
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the importance of a particular isospin-/ state depends on the nature of the interaction 
Hamiltonian. 

Of course, the Pauh exclusion principle can be considered as an isospin symmetric 
interaction, which forces the quantum states to be totaly anti symmetric in permutation of 
identical particles. In nuclear matter, basic constituents are protons and neutrons, which 
are subjected to Pauli exclusion principle. To appreciate it transparently, let us assume 
that the interactions can be turned off, namely, free gas approximation. Since there is no 
interaction, the reduced matrix elements for /-channel lose the /-dependence which means 
that there is no composition dependence of the energy: 

&''^%Np, Nn) = E^''^%N) (10) 

where E is the energy of the free-gas system without Fermi statistics. However, when we 
apply the Pauli principle to the system, the energy for the symmetric state {Np = Nn) should 
be lower than the state with protons or neutrons only. So the composition dependence is 
present even in the free Fermi gas system: 

Ef^^^iNp = N,Nn = 0)> Ef^^^iNp = 1/2, = 1/2). (11) 

We will return to this below. This aspect together with the isospin-depcndent nuclear force 
will be exploited for a transparent description of the symmetry energy. 

2.2 Nucleon-nucleon interactions 

For a system of N nucleons, the energy differences of the states with different compo- 
sition of protons and neutrons are encoded in the symmetry energy defined by subtracting 
the energy of the state with symmetric compositions, Np = Nn = N/2, from the energy of 
system composed of Np protons and Nn neutrons, 

EsymiNp, Nn) = E{Np, Nn) - E{Np = N/2, Nn = N/2) (12) 

or 

Esym{N, x) = E{N, x) - E{N, x = 1/2) (13) 
with X = Np/Nn- Clearly Esym{N,x = 1/2) = and 

Esym{Np, Nn) = Esym{Nn, Np), (14) 

hence 

Esym{N, X) = EsymiN, 1 - x) . (15) 

Expressed in terms of the energy density e = E/V, we have 

e(n, x) = e(n, x = 1/2) + esym{n, x), (16) 
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where n = N/V is the nucleon number density. Written in the quadratic approximation 
used in the hterature^ 



^sym 



{n,x)-^nS{n){l-2xf. (17) 



One can see this trivially holds for non-interacting systems, say, within 5 % accuracy, as 
given by 



47r^{n,x) = nSf^'^^{n){\-2x)\ (18) 



where 



5/r-ee^(22/3_l)£^0^iL)2/3 (19) 

5 no 

and^;0=_I.(^na)^/^ 

In what follows, we will focus on the "symmetry energy factor" (or simply symmetry 
energy) S{n) 

S{n) = [e(n,0) - e{n,x = l/2)]/n. (20) 

Here we assume that Eq. (17) holds^ for any value of x, for example see [16] [17]. 

The overall structure of the symmetry energy can be inferred by considering the two 
extreme cases with x = and x = 1/2. The configuration with x = with no protons is 
governed by neutron-neutron interactions only, and similarly for x = 1 with no neutrons but 
only with protons. Consider two identical systems of neutrons only with nucleon number 
N/2 (or density n/2) for each. Define the sum of the two e{n,x = 1/2) given by 

e(n, x = 1/2) = 2e(n/2, x = 0). (21) 

Now consider a nuclear matter consisting of N/2 neutrons and N/2 protons, for which the 
energy density is given by e{n,x = 1/2). Then A„p defined by 

Anpin) = e{n,x = 1/2) -i{n,x = 1/2) (22) 
= e{n,x = l/2)-2e{n/2,x = 0) (23) 

can be used as a measure of the n — p interactions. Then we have for the symmetry energy 
for the X = configuration with density n 

esym{n,0) = e(n,0)-e(n,x = 1/2), (24) 
= e(n,0)-2e(n/2,x = 0)- A„p(n). (25) 



^Higher order corrections - in (1 — 2x) - are known to be small near no but could be substantial at high 
density, ~ 7% of the quadratic term tsym at p > 6no[15]. 

^There may be deviation from this in the presence of kaon condensation. See below. 
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Defining A„„(n) as a measure of what could be called "double-neutron" interactions by 



A„„H = e(n,0) - 2e(n/2,x = 0) 



(26) 



we have 



^sym{n, 0) = Ann{n) - A„p(n) . 



(27) 



From eq.(17), we have for the symmetry-energy factor 



S{n) = [Ann{n) - Anp{n)]/n. 



(28) 



This is the expression wc will use for our discussions that follow. 

One can make a few simple observations with this formula. For the system of non- 
interacting nucleons, A„p(n) = 0, but A„„(n) > using the free Fermi-gas approximation. 
As is well known, this gives a positive symmetry energy factor S. Suppose that n-p inter- 
actions are attractive, i.e., Ajjp(7T,) <C 0. Then the "attrcLction" in the intevaction incveases 
the symmetry- energy factor, that is, makes it effectively "repulsive. " For those outside of 
the field, this aspect can cause confusion.^ 

3 Anatomy of the Symmetry Energy 

There are many approaches to the symmetry energy in the literature which can be 
broadly categorized by (a) phenomenological and (b) effective field theoretic (EFT). We 
shall now use Eq. (28) to unravel how nuclear forces control the symmetry energy. 

3.1 Phenomenological approaches 

As the first example, we consider a phenomenological model for the energy density 
and symmetry energy factor S.^ Eschewing going into details of pros and cons of various 
models, we take one simple approach called "MID" (short for "momentum-independent 
interaction") used by Li et al. [18] referred below to as LCK. It has the generic form that 
illustrates our basic point. Even though by definition, phenomenological models mix in 
basic interactions with many-body correlations, so it is difficult to single out symmetry 

For instance, at normal nuclear density TlQ , /Ann (n) includes the effects of "repulsion" due to Pauli 
principle. It gives ~ 13 MeV. The attractive n-p interaction, contributing with the opposite sign, gives the 
"repulsion" A„p(n) ^ 18 MeV, making the total ~ 30 MeV as is observed in nature. Roughly, A„„(n) is 
determined by the Fermi-Dirac statistics and A„j,(n) by the n-p interactions. More quantitatively, however, 
A„p(n) and A„„(n) depend strongly on the details of the nucleon-nucleon interactions as we will illustrate 
next. 

^The variety of sophisticated nuclear many-body approaches - including three-body forces - with pa- 
rameters adjusted to available experimental data, can describe the EoS up to the density constrained by 

experiments. However going beyond no involves extrapolations constrained neither by experiments nor by 
theory. The purely phenomenological model with adjustable parameters that is considered here encompasses 
such many-body approaches, so we will not discuss them here. 
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effects from dynamics, it is feasible to see how they enter into the symmetry energy. For 
this, we rewrite Eqs. (5.1) and (5.5) of [18] as 

e-*(n,x = l/2) = n[^- + B{-yi (29) 

S'^^'in) = Fa— + {18.6- F^){—f-, (30) 
no no 

where A,B,Ca,Fa,j are the parameters to be determined by experimental data. The 
energy density e and symmetry energy factor (or symmetry energy in short) S are now 
given by 

e{n,x) = n[mM + ^EU-?/' + ^- + Bi-r 
5 no 2 no no 

+{l-2xfS{n)], (31) 

Sin) = (2'/^-l)^EU-f/' + Fa^ + {18.6-Fa)i^f'^. (32) 
5 no no no 

The parameters A, B and 7 determined by experiments are A = —298.3 MeV, B = 110.9 

MeV, 7 = 1.21 and = 37.2MeV. The parameters F and C, however, cannot be uniquely 

determined. Thus there are possible values for F and C labeled by a (denoted x in [18]). 

Then using (22) and (26), we get 

A„„ = n[(2V3 _ i)|^o (iL)2/3 
o no 

4 no no no 

+(18.6-F„)(— )^«(l-2-^«)], (33) 
no 

1 T? 77 It? 

A„, = n[-^(-) + i3(-)^(l-2-^)--F„(-) 
4 no no / no 

-(18.6-Fa)(— )^°2-^"]. (34) 
no 

We can see the nature of degenerate Fermi gas only appears in A„„ as expected from the 
forms of the parametrization for the energy of the symmetric matter and the symmetry 
energy. 

The results for n = no for three values of a's, a = —1, 0, 1, are summarized in Table 1. 



a 


-1 





1 




-18.5 


-23.1 


-27.9 


^nn / n\no 


13.2. 


8.6 


3.8 



Table 1: A in MeV 



Recalling that Af^^^/n\nQ = 13.1 MeV for the free nucleon system without interac- 
tions, we see that while the attraction in the n-n attraction is negligible for a = —1, it 
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becomes much stronger for a = 1, with the exclusion repulsion largely compensated by the 
attraction. As for A„p, the n-p attraction becomes stronger, the greater the a parameter. 




Figure 1: The symmetry energy factor S{n) (thick line), A„p/n (dashed), A„„/n(dash- 
dotted) and the symmetry energy for non-interacting nucleons (dotted) for the LCK model 
with a = —1 (left panel), = (center panel) and = +1 (right panel) referred to as "super- 
soft." The vertical axis is in unit of MeV and the horizontal axis is n/no. 

The density dependence of A„„ and A„p are shown in Fig. 1 for a = — 1, 0, 1. These 
results illustrate how strong interactions via nuclear forces can influence the S - an indis- 
pensable part of the EoS, giving an answer posed in Introduction. While fine-tuned to fit 

experiments up to ~ no, the structure of a = 1 for n-n interaction is basically different 
beyond uq from both a = and a = —1. As we will discuss below, the key role here is 
played by the nuclear tensor force. What governs S above no has therefore a huge influence 
on the structure and the fate of compact stars. 

3.2 Relativistic mean-field approaches 

An alternative to the strict phenomcnological approach anchored on standard nuclear 
many-body approach is to resort to phenomcnological Lagrangians that incorporate as 
complete relevant degrees of freedom as feasible for practical calculations. One writes an 
effective Lagrangian in the spirit of doing effective field theory and treat it in the mean field 
ignoring higher-order loop corrections embodied in chiral perturbation theory, and fixing 
the parameters of the Lagrangian from experiments done in zero-density vacuum. We 
are therefore ignoring possible nonlocality associated with loop graphs. The Lagrangian 
is chosen on the basis of "naturalness" imposed on high order terms consistent with the 
power counting prescribed by the symmetry involved, chiefly chiral symmetry. The basic 
premise of such an approach is that in doing nuclear physics, the mean-field theory with 
an effective Lagrangian with correct symmetries is equivalent to doing Landau Fermi-liquid 
theory for nuclear matter [19] and that such a mean field theory can be extrapolated to 
high density regime relevant to compact stars. We shall call this "standard relativistic mean 
field approach" (SRMFA for short) to distinguish it from the mean-field of the hidden local 
symmetry (HLS for short) Lagrangian that we will discuss later. In the next subsection we 
shall propose that hidden local symmetric approaches implemented with both scalar and 
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baryonic degrees of freedom and appropriate scaling of the parameters is more predictive 
when treated in mean field. 

Here we will briefly illustrate what one can expect and what one has found in standard 
phcnonicnological mean-field approaches. In considering symmetric and asymmetric nuclear 
matter in mean field, we restrict to two flavors and ignore the pion field. The relevant degrees 
of freedom for the case at hand are the nucleons N, the p and co mesons and a low-lying 
scalar that we will denote s. It is important that the s is identified as a chiral scalar, not 
the fourth component of the chiral four vector usually denoted a. The leading Lagrangian 
with the lowest-dimension field operators take the usual form which can be read off directly 
from the hidden local symmetric Lagrangian (41) in Section 3.3.1 that contains both the 
dilaton scalar x ~ which wc will identify later with s - and the nucleon field N. What is 
not found in (41) but added in SRMF are higher-dimension field operators of the form 

SjC{N,S,lo^,P^) = bsiQssf + Caigss)"^ + c^iglu^oo'^f 

+9Xp'"'{dsgy + dvgloo.Lon + ■■■ (35) 

where gs,uj,p and b, c, d arc arbitrary constants fit to experiments. As an illustration, we 
take the parametric forms discussed by Piekarewicz and Centelles [49] up to second order 
in the density fiuctuation. The energy density of the asymmetric matter is written as 

^srmfin, x) = n[mN + eo + ^^ou^ + S{1- 2x)^] (36) 

where the symmetry energy factor is 

S = J + Lu+ ^KsymU^ . (37) 

Here the quantities J, L and Kgym are given in terms of the parameters of the Lagrangian 
and mean-field values of S, po S'ld ojq - the explicit forms of which we need not write 
down, u = ^^3^'^°^ , and uq = 0.148 fm~^^ corresponds to the equilibrium density. In 
this parametrization, the kinetic parts are not explicitly shown in contrast to LCK. The 
parameters, J,L,Ksym, the energy density of the symmetric matter cq and the compression 
modulus Ko are given in Table 2 for two models "FSU" and "NL3" both of which are fit to 
reproduce what is given by experiments up to ~ no. Prom (36) and (37), we find 

Ann = n[^Ko{u'^ -u^) + L{u-u) + ^Ksyrn{u'^ -u"^)] (38) 
Anp = n[^Ko{u'^ -u'^)-{J + Lu+^KsymU% (39) 

where u = ■ The resulting symmetry energy is shown in Fig. 2 for the parametric 

forms of FSU and NL3. 

Put together with the results of the standard phenomenological nuclear models, one 
can say that while constrained by nature up to ~ no, the symmetry energy for n > no 

®This is somewhat smaller than the canonical value 0.16 fm~*. 
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Model 


J 


L 


K 


Ko 


eo 


FSU 


32.59 


60.5 


-51.3 


230 


-16.30 


NL3 


37.29 


118.2 


100.9 


271.5 


-16.24 



Table 2: Bulk parameters characterizing the energy of symmetric nuclear 
matter and symmetry energy [49]. Kq and eg are, respectively, the com- 
pression modulus and the energy density of the symmetric matter at its 
equilibrium. 




Figure 2: Same as Fig. 1 for the FSU model (left panel) and the NL3 model (right panel). 



is totally undetermined. There is however one common feature in all these approach and 
it is that the attractive n-p interactions (A„p < 0) become repulsive (A^p > 0) at higher 
densities as seen in Figs. 1 and 2. 

3.3 Effective field theory approaches 

3.3.1 Hidden local symmetry 

Given the wildly diverging predictions by the available models, it is clearly desirable 
to devise theoretical controls on the relevant quantities. At present no model- independent 
field theoretic approach is available for dense baryonic matter^: lattice QCD cannot access 
that regime. The approach closest in spirit to - though not directly derived from - QCD 
is the skyrmion description of dense baryonic matter, and the first attempt to deduce the 
symmetry energy from it was made in [21]. Anchored on large Nc properties of QCD, 
the approach captures nonperturbative aspects of large Nc gauge structure of the strong 
interactions - and most importantly for any number of A baryons for A = 1, oo with A 
identified with the topological winding number B of the soliton. While the skyrmon picture 

'^There are attempts to apply chiral perturbation theory (BTrChPT) with baryonic chiral Lagrangian 
(involving nucleons and pions only) to nuclei and nuclear matter. BTrChpT has some modest success in 
nuclei and nuclear matter but does not work for some quantities that are intricately tied to Fermi-liquid 
fixed-point structure of nuclear matter such as the orbital gyromagnetic ratio Sgi for which mean-field 
approach works better [20]. In particular, it is not likely to work for dense matter we are concerned with. 
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was introduced by Skyrme for baryons [22], the concept has proven to be extremely powerful 
and pervasive in many areas of physics, most successfully in condensed matter physics [10]. 

In applying this concept to dense baryonic matter, however, there is a technical 
difficulty and that is, it is not known how to write down an effective Lagrangian that 
captures QCD in a manageable form. Written in the most minimal way in the large Nc 
limit is the Skyrme Lagrangian that contains pions only, consisting of the two-derivative 
current algebra term and the quartic term (called Skyrme term) that provides topological 
stability to the skyrmion. There is nothing in principle that prevents from bringing in 
higher derivative terms as the energy scale is increased but one cannot work with such a 
Lagrangian that involves too many unknown parameters to control. 

An alternative field theoretic approach is density-functional approach anchored on 
mean-field treatments of effective Lagrangians with pertinent symmetries taken into consid- 
eration. This approach is essentially related to the standard relativistic mean-field (SRMF) 
theory mentioned in Section 3.2. Here again the plethora of higher-order terms involved in 
the approach are difficult to control systematically. 

One way to circumvent this difficulty in both approaches is to bring in vector de- 
grees of freedom that figure at higher mass scales. This can be done most effectively by 
exploiting hidden local symmetry which can appear in the leading Nc order on the same 
footing with the pion. This can be done naturally by exploiting the redundancies in- 
herent in the chiral field^ U = e^^'^l^- G SU{Nf)L x SU{Nf)R/SU{Nf)v=L+R, namely 
^ = ^L^R = ^L9H^)9i^)^R where one has g{x) G SU{Nf)v- Gauging the redundancy g{x) 
by introducing the gauge fields V makes the nonlinear sigma model Lagrangian extended 
to the energy scale of V. The gauge fields here can be considered as "emergent" from pion 
interactions going up from the current algebra scale. For Nf = 2, the lowest-lying hidden 
(or "emergent") gauge mesons are p and uj, so we will be concerned with V = p,u;. 

The power of hidden local symmetry is that it allows a systematic chiral counting, 
hence a systematic chiral perturbation expansion with vector mesons included. Further- 
more it can consistently treat the situation where the vector meson masses become as light 
as the pion mass - which is massless in the chiral limit, a feature which is lacking in phe- 
nomenological models with no local gauge symmetry [23]. 

Another power in this approach is that an infinite tower of vector mesons can be 
incorporated as hidden gauge fields. This can be derived bottom-up by exploiting the 
infinite redundancies and generating the infinite tower of hidden gauge fields, yielding a 
de-constructed 5D Yang-Mills action [24]. It can also be obtained top-down by reduction 
via branes from string theory. Dual - and closest - to QCD with correct chiral symmetry 
is the Sakai-Sugimoto hQCD model [25]. 

What the impact of the higher tower of vector mesons on RMF approach is has not 
yet been investigated. On the other hand, for the skyrmion approach, it is found that the 

*In this paper we will confine to two flavor SU{2). 
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tower of HLS fields encoded in the 5D YM action markedly improve on the structure of 
skyrmion for single baryon as well as for many-baryon systems [26]. By taking a flat space 
5D YM action and using the BPS skyrmion, Sutcliffe shows that the more gauge fields 
incorporated, the better the skyrmion description of the baryonic systems. This implies 
that broken scale symmetry tends to be restored as more vector mesons are incorporated 
into the skyrmion structure. We will see that this feature emerges in the skyrmion crystal 
we consider in this paper. 



3.3.2 HLS with baryons and dilaton 

In describing nuclei and nuclear matter, a low-lying scalar meson degree of freedom, 
either as a local field or composites, is essential. It gives the necessary intermediate-range 
attraction needed for the binding of the system. It is however missing in HLS with or without 
the infinite tower. How to incorporate scalar fields in hidden local symmetric Lagrangian 
is problematic because of the unnaturalness associated with scalar fields and has not been 
worked out satisfactorily in that gauge theory. Focusing on mean-field approaches, we 
choose to introduce the needed scalar degree of freedom via the trace anomaly of QCD 
and introduce the scalar as the dilaton that accounts for the "spontaneously broken" scale 
invariance in the presence of the explicit breaking due to the anomaly. How this can be done 
is explained in [27]. The key assumption in [27] is that there is a "soft" dilaton Xs that is 
locked to chiral symmetry the vev of which vanishes when chiral symmetry is restored and a 
"hard" dilaton Xh whose vev remains non-vanishing above the chiral critical point, playing 
the role of an explicit symmetry breaking to scale symmetry in the low-energy (confinement) 
sector as required to make the spontaneous symmetry breaking viable [28]. We identify Xs 
as the pseudo-Goldstone scalar in the spontaneous breaking of scale symmetry. 

One can readily write down such a dilaton-implemented HLS (dHLS for short) La- 
grangian [12]. Retaining only the lowest vector mesons p and u, and integrating out the 
"hard" dilaton, the dHLS has the form^ 

CdHLS = Cm + C^ , (40) 

Cm = i^X^ tr [a^j^a^] + aKX^ tr - ^tr [V^uV"] + 

1 KUl'i 

where we write x for x^, a = {F^/Ft^Y = ^^/^x^ "^x mass of the soft 

dilaton. The quantities and aff are defined in the review [23] and in [12]. The ellipsis 




Considering Nf = 2, we write the Lagrangian in the U{2) symmetric way with the vector field defined 
as gVi_i = g{^T ■ + ^t^J^t). The U{2) symmetry seems to work well in the zero-density vacuum, but there 
is no good reason for it to be a good symmetry in medium. We will however assume ^(2) symmetry up 
to the density at which skyrmions make the transition to half-skyrmions. For the symmetry SU{2) x (7(1), 
gVij, should be replaced by {^t ■ + ^ui/j,). 
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stands for higher order terms that result by integrating out the higher members of the 
tower. 

We will consider baryons arising from this Lagrangian or rather a truncated version of 
it as solitons in what follows, that is, in putting skyrmions on crystals. But it is interesting 
to put baryons explicitly into the Lagrangian as one does in baryon chiral perturbation 
theory. At a given order, the two ways of treating baryons, as skyrmions or local baryon 
fields, is equivalent and this equivalence was used also for describing intanton bayons in the 
SS' hQCD model [29, 30]. With both the dilaton field and the baryon (nucleon in this case) 
field, the Lagrangian has the form [12] 

^bdHLS = + + i^x ' i'^^) 

with 

Cn = NipN - ^mNNNx + gAN^±l5N + gvN^N + ■■■ . (42) 

Cm and are as defined above. This is for flavor U{2). For SU{2) x ^7(1), gyN^^^N 
should be replaced by {gvpN^^^N — gvoj^^N). 

3.3.3 bdHLS and the dilaton-limit fixed point 

We will discuss mean-field properties of the Lagrangian (41) below. Here we consider 
how it can describe both low-density matter - near the density of nuclear matter - and 
high density matter near chiral restoration. For suitably adjusted parameters such as the 
dilaton and vector-meson masses fit to data, we can employ the Lagrangian (41) in mean 
field to give nuclear matter. This is essentially Walecka linear mean-field model. With 
the parameters endowed with appropriate density dependence [40] , it can actually describe 
fairly well all nuclear matter properties [31]. The reason for this is that the mean-field 
theory so formulated is equivalent to Landau Fermi-liquid fixed point theory and therefore 
it should work well very near the nuclear saturation point. However there is no reason to 
expect that this theory extended in a naive way to high density would work. In fact near 
the chiral transition, the scalar dilaton which is a chiral scalar cannot be the relevant scalar. 
It must instead be the fourth component a of the chiral four vector {7^1,^2,1^3, a). How the 
chiral scalar can be transmuted to a by density is unknown and needs to be clarified. 

In [12], the solution was found in terms of what is called "dilaton limit" by Bcanc 
and van Kolck [32]. In what follows in this subsection, we will confine to the flavor SU{2) 
sector of vector mesons, i.e., p mesons, returning to the U{1) sector later^*^. Combining the 
dilaton field x and the chiral field U, redefine the fields as 

E = UxVi^ = £.URxVK = s + iT-TT, (43) 

^"Considering U{2) symmetry in [12] is not justified in medium since there is no reason to expect that the 
symmetry holds in medium as it does in free space. The discussion in [12] on the suppression of repulsion 
due to w-exchanges therefore may be incorrect. 
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+ +75(d-d)]^> (44) 



with the Pauh matrices r in the isospin spaee. Expressed in these fields, the "linearized" 



Lagrangian includes terms which generate singularities, negative powers of tr 



that 



goes to zero in chiral symmetric phase. Those terms appear as multiplicative factors: 

Xn = gvp- 9A, Xy^ = l- K. (45) 

Here gvp represents the "induced" /O-N coupling replacing the gv relevant to C/(2) sym- 
metry in [12]. In order to prevent the divergences, one is then required to impose that as 
tr [sStj ^ 0, 

Xn ^0, X^^ 0. (46) 

We assume that the large Nc approximation is reliable for gA at large density which 
would imply that — >^ 1 as the density increases toward chiral transition. At the point 

K = gA = gvp = ^ (47) 

the bdHLS Lagrangian goes over to Gell-Mann-Levy sigma model. This is referred to as 
"dilaton limit" (DL for short). The direct vector- meson coupling to the baryon in (41) 

9VN = 5p(l - gvp) (48) 

where gp is the SU{2) hidden gauge coupling therefore goes to zero 

gvN ^ 0. (49) 

Thus approaching the DL, the vector mesons decouple at the tree level from baryons al- 
though V-TT couplings remain. 

An interesting possibility with this limit is that it can be a fixed point of the bdHLS 
theory. Let us call it dilaton-limit fixed point (DLFP for short) and the density at which 
the fixed point is reached will be denoted ridi- Indeed it has been shown that DL is an 
immediately visible fixed point in the one-loop RGE of (40) [13], This theory may have a 
variety of other fixed points as in the original HLS theory [23] out of which the vector mani- 
festation fixed is the only fixed point matched to QCD. Now we conjecture that approaching 
the Gell-Mann-Levy sigma model is tantamount to approaching the chiral restoration point, 
and that the DL is the fixed point of QCD in the same spirit as the VM fixed point is. It 
remains to show that there are no other fixed points in this bdHLS theory (41) that can be 
matched to QCD, in the same vein as that there are no other fixed points than the VM in 
HLS that match to QCD [23]. 



^^This is seen in lattice calculations. It is found that qa calculated in quenched approximation in lattice 
QCD - which corresponds to the large Nc limit - is close to what is found in unquenched calculations. 
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One interesting question that remains unanswered is whether the DLFP and the VM 
fixed point in HLS theory are at the same density or even related to each other. It has 
been established that the VM fixed point remains intact in the presence of fermions whose 
dynamical mass goes to zero (in the chiral limit) at the chiral restoration point [33]. This 
would be the case if baryons were chiral with masses generated dynamically. It seems highly 
likely that it will be the same if the baryon can be treated as heavy near the chiral transition 
point - as in the mirror assignment for the baryons with a large chiral invariant mass [12] - 
in which case the fermion will decouple when approaching the VM point. It is not obvious, 
however, what would happen if the fermion were neither light nor heavy. Although much 
less studied, the situation seems quite similar for the DLFP (47) as shown in [13]. As we 
will see, our analysis suggests that they cannot be too far apart. 

3.3.4 Half-skyrmion matter and BLPR scaling 

We now turn to the skyrmion description of dense matter - which is our principal 
tool for our problem. A work is presently in progress [34] but no results are yet available 
on skyrmions in dense matter anchored on the dilaton-HLS Lagrangian advocated above. 
We shall therefore pick a vastly simplified model and extract relevant information from it 
in applying to nuclear and dense matter. We will not be able to obtain numerically precise 
results, but we consider them to be qualitatively robust and then develop a procedure to 
make quantitative calculations. 

The model we shall use is the Skyrme model that includes the soft dilaton field. We 
shall put the skyrmions of this model on the crystal and deduce the mean field features of 
the ground state of baryonic matter. 

We first recall what one obtains in the absence of the dilaton, that is, the pure Skyrme 
model. From the skyrmions put on an FCC crystal with the chiral field described by Atiyah- 
Manton's holonomy ansatz [35], we learn [36] that in the half-skyrmion phase to which 
skyrmions make transition at a density 7 scale invariance, broken in the skyrmion phase, 
gets restored. This can perhaps be understood as resulting from a "symmetry enhancement" 
associated with lattice translation [37]. Here scale symmetry and chiral symmetry are 
connected in such a way that the restoration of scale symmetry is locked to the vanishing 
of the chiral condensate {qq)- In this model, the half-skyrmion matter is in the phase where 
the quark condensate {qq) vanishes on the average in the unit shell, but locally non-zero, 
thereby producing an inhomogeneous structure. 

The situation is different when the dilaton is present. Consider adding the soft dilaton 
to the Skyrme model of [36] that accounts explicitly for spontaneously broken scale symme- 
try which gets restored at the chiral transition density n^. We shall call this "Dskyrmion" 
model. The Dskyrmion Lagrangian is constructed such that the scale symmetry is sponta- 
neously broken by a Coleman- Weinberg-type potential V{'x)- It can be obtained from the 
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dHLS Lagrangian (40) by dropping the vector mesons and inserting the Skyrme term, 
Csk = ^(l) Tr{L,Ln + ^Tr[L„L^f 

+^ (j-^ TiM{U + U^ -2), 

+ld^xd^X + V{x) (50) 

where Ai is the quark-mass matrix, V{x) is the potential that encodes the trace anomaly 
involving the soft dilaton, = 11^8^11, with U the chiral field taking values in su{2) and 
is the vev of x- 

The results of the analysis [21, 38] are as follows: 

1. The skyrmion-half-skyrmion phase transition takes place at a density ni^2 slightly 
above uq. The precise value of ni/2 does not depend sensitively on the dilaton mass 
for < 1 GeV. For a light dilaton mass appropriate for the soft component, i.e., 
~ (600 — 700) MeV, which is also reasonable for nuclear phenomenology in the vicinity 
of nuclear matter density, the transition takes place in the range 1.3 < ni/2/'^o ^ 2. 

2. At ni/2, the quark condensate {qq) oc Trll — ^ but the pion decay constant F* 
remains non-zero, dropping slowly at increasing density up to the chiral transition 
density riy^, signalling that the half-skyrmion phase is in the Nambu-Goldstone mode. 
The chiral restoration density is on the contrary highly sensitive to the value of 
the dilaton mass. For the mass adopted here, it comes out to be > 4no. 

3. In the half-skyrmion phase, the nucleon mass scales with the dilaton condensate xo = 
ix) ~ which drops slowly, not with the quark condensate. This is similar to the 
behavior of the technibaryon discussed recently as a source for dark matter [39] . It 
remains non-zero up to the critical density n^. The dilaton mass scales also linearly 
in (x). 

The phase structure given by the Dskyrmion model is then as follows. For (a) < 
n < ni/2, the system is in the skyrmion phase, i.e., in the standard Nambu-Goldstone mode 
with a non-vanishing quark condensate and pion decay constant; for (b) ni/2 < n < n-^, 
the system is in the half-skyrmion phase with the vanishing {qq) but non-zero pion decay 
constant. Chiral symmetry is still broken with an order parameter not characterized by 
{qq) but by a vev of higher dimension operators; for (c) n > n^., the system is in standard 
Wigner-Weyl phase with vanishing quark condensate and pion decay constant. Our model 
cannot access this phase. 

This half-skyrmion structure at ni/2 is expected to bring about substantial changes 
at densities above normal nuclear matter density in hadronic properties, in particular, the 
Brown-Rho (BR) scaling proposed in [40] . In this paper we focus on nuclear tensor force in 
which scaling has a strong influence as first argued in [41]. 
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In order to write the new scaling (dubbed "BLPR" for short), we need to bring in 
the vector degrees of freedom. The pertinent Lagrangian is the dHLS (40). Up to date 
no consistent and trustful calculations putting this Lagrangian on crystals arc available . 
We can however exploit that HLS is gauge-equivalent to the non-linear sigma model and 
deduce their scaling based on what is known of HLS. HLS predicts that approaching the 
chiral transition at n^^, the vector meson mass scales as [23]^^ 

m*/mp ~ 5p/5p ~ (^^')7 (w)- (51) 

This implies that the p mass vanishes as the quark condensate goes to zero. If the local field 
description for meson degrees of freedom makes sense as one approaches chiral restoration, 
this would imply vanishing mass at the chiral phase transition density n^. 

If the flavor U (2) is not a good symmetry in medium - which is likely the case at a 
density n > ni/2 ~ for the p and the scaling of the lo mass need not be the same as that 
of the p. The vector manifestation may not apply to the singlet vector meson. However it 
seems plausible on the ground of "mended symmetries" that the u mass also goes to zero. 
How it does is not known. 

It has been shown that (51) holds in the presence of baryons whose masses are entirely 
dynamically generated [33] (with the baryon mass going to zero at n^). If on the other hand, 
the baryon mass remains massive as was found in the Dskyrmion model discussed above 
or the parity-doublet model with a large mo [12], the baryons can be integrated out in 
approaching n^, so we expect that the scaling (51) should still hold. 

The PLBR we shall study is the following. 

• For the "low density regime (regime-I)" < n < nxji , the scaling remains the same 
as in the old BR [40]: 

m*j^/mN ~ niy/mv ^ m*/ms ^ F*/Ft, = ^i{n), 
gv/gv - 1 (52) 

^^A systematic work is in progress on this problem. There have been performed several calculations 
(involving one of the authors (MR)) with HLS Lagrangian with or without the dilaton scalar. Unfortunately, 
however, none of them is fully consistent with the symmetries involved: The homogeneous Wcss-Zumino 
terms that are present in the presence of vector mesons are not correctly treated. We therefore do not rely 
on the results so obtained. 

^^We have not verified that one can approach the vector manifestation (VM) fixed point at which chiral 
symmetry is restored from the half-skyrmion phase. Hence the scaling (51) has not been proven on crystal 
lattice. We are simply assuming it. One can however surmise what could happen. As stated above, the quark 
condensate (qq) is zero in the half-skyrmion phase while with the dilaton field present, chiral symmetry is 
not necessarily restored. One way to interpret the scaling (51) in the context of skyrmion crystals is that 
in the half-skyrmion phase, (qq) ~ TiU = TtUi + TtU2 ~ when averaged over the cell with the two 
half-skyrmions giving Tr[/i = c and TrC/2 — —c with c 7^ and that as one approaches the VM fixed point, 
c 0. Possible non-vanishing chiral order parameters on the crystal - conjectured above - could carry this 
quantity c in such a way that they vanish at n^. 

^''We have replaced the dilaton x by the scalar s that represents the scalar that figures in nuclear physics. 
It will be mostly of x at low density but could be different at high density with possible admixing to glueballs. 
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where V stands for both p and uj. 

• For the "high density regime (regime-II)" 77,1/2 < n < n^, the scahng is changed to 

mlf/mN ~ b, 
m*/mp « g*/gp = ^ii{n), 

ml/ms ~ ^s{n). (53) 

Since the nucleon mass is found to scale weakly, we set h equal to a constant h<l. ??? 
There are four scaling functions ^iji^uj,s that can be different though they must be 
related to each other. The bdHLS model supplemented with the dilaton-limit and VM 
fixed points should provide their behavior in density. For instance, although both the 
nucleon mass and the dilaton mass, at first sight, scale with the dilaton condensate 
(x), we can see from the analysis in [42] that the scalar mass drops much faster 
than the dilaton condensate at density > hq. There are several reasons to expect that 
<I>s(n) < b. First there seems to be a strong back-reaction from the soliton background 
on the fluctuation of the dilaton field whereas there can be a substantial-sized chiral- 
invariant mass mo in the nucleon mass as in the parity-doublet model [12, 13] which 
prevents a rapid drop. More generally, locked to chiral symmetry, must join the 
pion mass m^r = at in the chiral limit. 

As for the scaling if U{2) symmetry held for p and cu, we would have = $77. 
We expect it not to hold in medium for n ^ ni/2 and hence we have no guidance from 
the vector manifestation, the only possibility that we have being that it vanishes at 
the chral restoration density.. 

Given that we have no reliable numerical results for the scaling functions, for numerical 
analyses, we shall simply choose reasonable values for Cx in the form 

for X = I, II, s, oj. 

3.3.5 Tensor forces 

A simple prediction as to how tensor forces behave in nuclear and dense matter can 
be deduced from the above scaling relations [21]. Since the nucleon mass remains more or 
less unsealed, we can use non-relativistic approximation for the nucleon and write down the 
tensor forces contributed by the tt and p exchanges. The well-known formulas are 

^m(0 = SM^p^mMTl ■ T2S12 



47r 

1 1 

+ 7 ^ + 



{■niMr)^ {niMrY 'iuiMV, 



(55) 
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where M = vr, p, Sp(^) = 1). Note that the tensor forces come with an opposite sign 
between the pion and p tensors. 

Following the common practice, we will leave the pion tensor unaffected by the 
density. ""^^ As for the p tensor, apart from the scaling mass m*, it is the scaling of /jvp 
that is crucial. Written in terms of the parameters of the bdHLS Lagrangian (41), we have 
the ratio 

^ fNp ^ gh^^m^ (56) 
fNp gvN rup m% 

where we have defined the effective V-N coupling gvN = 5(1 ~ 9v)- It follows from the 
scaling relations (52) and (53) that 

« 1 for 0<n<ni/2 (57) 



$7/ for ni/2 < n < nc- (58) 



We note that the ratio R will be strongly suppressed for n > ni/2- In addition to the 
suppression by the factor the approach to the dilaton-limit fixed point would bring 
(1 — gy^ go to zero. This makes the p tensor killed rapidly. 

What happens is illustrated in Fig. 3 both with the "old BR" (no topology change) 
and the "new BR" for BLPR (with topology change). With the old BR where the skyrmion- 
half-skyrmion phase change is not taken into account, the net tensor force decreases contin- 
uously in density with the attraction vanishing at n 3no for the given parameters. In a 
stark contrast, with the BLPR, the p tensor force nearly disappears at n ~ 3no leaving only 
the pion tensor active. Up to nij2, they behave the same way, so that the C14 dating beta 
decay [43] will be unaffected by the BLPR. Once the density exceeds ni/2) there is a drastic 
difference between the two scenarios. In the BLPR, the p tensor gets rapidly suppressed 
and for a reasonable values of c^, the pion tensor is the only component left. This will have 
an interesting consequence on the nuclear symmetry energy as mentioned below. 

3.3.6 Effect of BLPR on the symmetry energy 

If we assume as argued by several authors [44, 45, 46] that the symmetry energy is 
dominated by the tensor force, then the main contribution to the symmetry energy factor, 
following Brown and Machleidt [44] , can be written as 

S - {Vsym) ^^^{V:^{r)) (59) 

where E ^ 200 MeV is the average energy typical of the tensor force excitation and Vr 
is the radial part - with the scaling factor R taken into account - of the net tensor force. 



^®In fact it can be shown with the help of low-energy theorems that the pion tensor force remains little 
modified up to n ~ (3 — 4)no. 
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Figure 3: Sum of tt and p tensor forces = (ti • T2Si2)~^{y^ + Vj) in units of MeV for 
densities n/uQ =1, 2 and 3 with the "old scahng" — 0.15n/no and i? Ri 1 for ah n 

(left panel) and with the "new scaling," k. 1 — 0.15n/no with R k. \ ion n < Uxj^ and 
$7 and R ^ for n> ni/2, assuming no < ni/2 < 2no (right panel). 



There are other contributions such as the kinetic energy term etc. They can be ignored for 
our purpose. We see immediately from the above discussion that the symmetry energy will 
decrease toward ni/2 and then increase afterwards as the cancelation diminishes, giving a 
cusp at ni/2- As described below, this cusp feature is supported by a pure skyrmion neutron 
matter when collective-quantized. 

This prediction is in contrast with that of the old BR discussed in [47] . In this article, 
the authors apply the old BR to explain the "supersoft" symmetry energy which vanishes 
near n ^ Suq, namely, the a = +1 curve in Fig. 1. This comes about because the net 
tensor force decreases beyond n ~ no in the old BR. In this paper we are arguing that the 
decrease is stopped by the presence of the half-skyrmion phase. Experiments will decide 
which scenario is viable. 
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3.3.7 Symmetry energy in the skyrmion-half-skyrmion phase transition 

In this subsection, we discuss the eym obtained [21] by quantizing the skyrmion mat- 
ter described by the Lagrangian (40). This confirms the correctness in transcribing the 
skyrmion crystal result into an effective Lagrangian. 

In the skyrmion framework, the symmetry energy comes from a term subleading in 
Nc- It must therefore arise from the collective quantization of multi-skyrmion systems. In 
his original work on the skyrmion crystal, Klebanov [48] discussed how to collective-quantize 
the pure neutron system. We apply this method to the skyrmion matter as well as to the 
half-skyrmion matter we have obtained. 

Consider an ^-nucleon system for B = ^ — ^ oo (where B is the winding number). 
Following Klebanov, the whole matter is rotated through a single set of collective coordinates 
U{f,t) = A{t)Uo{f)A^{t) where Uo{r) is the static crystal configuration with the lowest 
energy for a given density. The canonical quantization leads to 

E^ot = AM^i + (^*°* + 1) , (60) 

where M^j and A/ are, respectively, the mass and the isospin moment of inertia 

which is given as an integral over the single shell of a function consisting of the 

skyrmion configuration Uq- 

Since is the total isospin which would be the same as the third component of 

the isospin I3 for pure neutron matter, this suggests taking^^ iov S = {N — P)/{N + P) < 1 

jtot ^ (gj) 

Thus the energy per nucleon in an infinite matter (A = 00) is 

E = Eo + -^^S\ (62) 

with £^0 = Mci- This leads to the symmetry energy factor 

S«±-. (03) 

Again the kinetic energy contribution (which is l/7Vc-suppressed and expected to be unim- 
portant at high density) is ignored in this expression. The numerical results for the given 
parameters are plotted for densities below and above no in Fig. 4. What is noticeable is 
the cusp at which is similar to what is expected with the BLPR as discussed above. 
This feature seems to be considerably smoothed out by nuclear many-body correlations in 
microscopic nuclear field theory approach with Viowk and old BR scaling incorporated into 
renormalization-group fiow. It leaves, however, its distinctive imprint in the change in the 
symmetry energy from soft to hard at the density of topology change. 



^^This is equivalent to assuming that Ci in Clebsch-Gordan series, Eq. (4), is dominated by 7 = 1/2|A^„ 
Np\. It is definitely the case with iV„ » Np ot Np » Nn- 
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Figure 4: Symmetry energy given by the collective rotation of the skyrmion matter with 
/tt = 93 MeV, 1/e^ 0.03 and two values of dilaton mass. The cusp is located at ni/2- The 
low density part that cannot be located precisely is not shown as the collective quantization 
method used is not applicable in that region. 

A few caveats with this prediction are in order here. 

First of all the Klebanov collective quantization method is not expected to be appli- 
cable for very low densities, so one cannot take seriously the results below the density no- 
Given that the Dskyrmion model (that is, the Skyrme model with dilaton) for the energy 
eo(n) at the classical level cannot describe nuclear matter saturation, one may ask how one 
can trust the calculation of egym- The possible answer to this question is that what we are 
calculating for egym is a 1/Nc effect like the N-A mass difference since it is oc 1/A/ and the 
moment of inertia Xj on Nc- What enters into A/ is the leading Nc term and the subleading 
effects that figure in the saturation of symmetric nuclear matter given by the energy Eq 
would not affect egym to the leading order we are considering. An evidence for this is the 
fact that the location of the density is extremely insensitive to the dilaton mass as one 
can see in Fig. 3 of [21] whereas the Eq - hence the density at which (^)* = /* = - is 
strongly affected by the dilaton mass. 

3.3.8 Dilaton-limit fixed point and mean-fi,eld theory 

The relativistic mean-field methods of the type mentioned in Sect. 3.2 assume that 
the parameters of the Lagrangian fixed in the matter-free vacuum remain unmodified when 
extrapolated to high density in the mean field. As mentioned above, doing mean field 
calculation with an effective Lagrangian endowed with pertinent symmetries for nuclear 
matter is justified because it is equivalent to Landau Fermi-liquid fixed point theory [20]. 
It is however not clear that it applies equally (1) to higher density departing far from 
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the Fermi-liquid fixed point and to (2) the number of flavors greater than two involving 
strangeness in addition to up and down quark flavors. 

Now in the standard RMF approach, possible modifications of the Lagrangian param- 
eters due to the presence of medium are incorporated by the mean-field of higher-dimension 
operators. In the "sliding-vacuum" adopted in this paper, such high-dimension effects are 
albeit approximately captured in the density dependence of the parameters of the effective 
Lagrangian with hidden local symmetry [20, 31]. Suppose one calculates the symmetry 
energy using then the Lagrangian (41) in the mean field. We get 

5/» = g+- (64) 

where gpM = ^(1 — gvp)- Here the ellipsis contains other contributions such as the kinetic 
energy term that we may ignore at high density. We think that the effective p mass m* 
encapsulates most, if not all, of corrections from higher- dimension operators in the standard 
RMF Lagrangian. 

Now let us see what happens when the system approaches n^. The p mass goes like 
~ g*p. Therefore S will go like 

S/n^{l-g*ypf^Q (65) 

due to the dilaton limit fixed point (49). This highly non-perturbative feature encoded in 
the DLFP is absent in the standard RMF approach or in finite-order chiral perturbation 
theory. 

This of course does not mean that the standard RMF is wrong in the relevant density 
regime. It would remain valid until the DL sets in, but that may be at higher density 
than relevant for the process. This could be the case if the DL were to set in at about the 
same density as the vector manifestation. What the DLFP means is that the symmetry 
energy must diminish at some high density within the hadronic phase. What happens at 
the cross-over (or phase-transitioncd) to a quark-gluon phase is unknown. It stops being 
relevant in any case once the quark-gluon degrees of freedom set in. 



4 BLPR and Compact Stars 

In this section we apply the formalism developed above to compact stars, addressing the 
issue of the maximum mass vs. radius of neutrons stars. This will be a concise summary 
of the work in progress reported in [50], to which we refer for references. 

We are interested in implementing the BLPR in an effective field theory (EFT) for 
nuclear matter. In doing this we adopt the "double-decimation" strategy [51]. 

1. The first decimation consists of obtaining via RG equation the V/o«)jfc in free space by 
decimating to the scale ^lawk that describes nucleon-nucleon interactions up to lab 
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momentum to ~ 300 MeV. It is best for our purpose to think of doing this in terms of 
a generahzed HLS Lagrangian with parameters with the intrinsic density dependence 
given by the BLPR. 

2. The second decimation is to do nuclear many-body calculation with this Vi^nuK to the 
Fermi-momentum scale Afermi- There are a variety of ways of doing this step. They 

all amount essentially to doing Landau Fermi-liquid fixed point theory and arrive at 
nuclear matter at equilibrium density no ^ 0.16 fm^. This is an EFT well-justified 
up to density near no- This step fixes the scaling properties of the parameters in the 
Lagrangian up to near n = uq. We extend the same scaling up to ni/2 since it is not 
too high above no- 

3. The last step in our application is to smoothly extrapolate with the formalism to high 
densities and calculate the EoS for compact stars. In doing this, one can adopt chiral 
perturbation strategy and include n-body forces - if one wishes - for n > 2, suitably 
introducing scaling parameters for the n-body forces. In [50], the topology change is 
incorporated in a smooth manner in terms of the changes in intrinsic scaling at ni/2- 
It however ignores other degrees of freedom that might enter, such as strangeness that 
manifests in terms of kaon condensation and strange quarks etc. 

What came out of this exercise can be summarized as follows. 

With the scaling properties constrained by nuclear matter properties in the region- 
I and heavy-ion data available up to n ~ 4.5no in the region-II, a topology change at 
1.5 < n/riQ < 2 can account for large- mass compact stars with no intervention of other 
degrees of freedom than nucleonic. A pure neutron-matter calculation gives M^ax ~ 2.4M0 
and RKi 11 km. The interior density of the star comes to ~ 5no. This gives a rough idea 
since compact stars are not pure neutron matter. A calculation that takes into account 
fully the relevant star conditions is in progress and will be reported in a later publication. 
Given that the density reaches the regime where kaon condensation can take place, possibly 
precipitated by the topology change discussed below, the appearance of strange quarks will 
also have to be addressed. 

5 New Degrees of Freedom: Strangeness and Quarks 

So far we have been concerned with the symmetry energy in the two-flavor sector 
with the relevant baryonic degrees of freedom given by proton and neutron. In compact 
stars, apart from leptons, other hadronic degrees of freedom are expected to enter. In this 
section, we will consider how kaons and hyperons can figure in the symmetry energy relevant 
to compact stars. For this, we relate the chemical potentials of the participating degrees of 
freedom to the symmetry energy. 
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5.1 Electron and muon threshold 



We first consider the proton-neutron system without strangeness in chemical equihb- 
rium in weak interaction with electrons (e) and muons (/x) . Denote the proton and neutron 
number densities by Up and n„, respectively. In terms of the chemical potentials 

dUr, dUp ^^^^ 

we have from Eqs. (16) and (17) 

IJ-n- IJ'p = 4(1 - 2np/n)S{n) = He = IJ'^ = 1^- (67) 
The system must be charge neutral, so 

np = ne + n^ = ^[^^ + e(/x - m^)(/x2 - m^)3/2] (68) 
From Eqs (67) and (68), we get 

This relates the chemical potential /i to the symmetry energy factor S for a given density. 
We will see how this chemical potential figures in the EoS for compact stars. 

5.2 Strangeness degree of freedom 

The first "exotic" degree of freedom that we consider is the anti-kaon, in particular 
K', that carries the strangeness. The pion could also figure but it is not likely to enter 
significantly into the process. We will ignore it here. Although in free space, (that we 
will simply refer to as "kaon" unless otherwise noted) is massive, its mass is expected to drop 
in dense medium because of strong attractive interactions with the nucleons in medium and 
when it becomes light enough, then the strangeness changing processes n p + K~ and 
e~ O K" + Ve can take place. When the kaon mass drops sufficiently low so that it crosses 
the increasing electron chemical potential jie at increasing density, kaons will condense and 
soften the EoS [52, 53, 54]. If these processes take place in chemical equilibrium with 
protons, neutrons, electrons and kaons in stellar matter, then 

fJ-n - IJ'p = fJ'e = fJ'K = IJ- (70) 

There have been a large number of discussions on the possible impact of kaon conden- 
sation in compact stars and we will not enter into details of the theoretical tools employed. 
See for a review [55]. We will focus here on certain novel aspects of the process that takes 
place in the hadronic sector that we consider to be of relevance but have not been so far 
discussed. We will come to the quark-gluon sector later. We shall also postpone to later 
the role of strangeness that manifests through hyperons. 
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To illustrate the basic point, wc sketch in the simplest form what is involved. For 
this we first consider chiral perturbation theory. Write the grand potential, SI of nuclear 
matter with s-wave kaons as 

n{n, X, n, K'^)/V = n{n, x, 0) /V - a{n, x,h)K'^ + --- (71) 

where we have expanded the grand potential up to quadratic order in the vev of the kaon 
field, the ellipsis standing for higher orders in K. To Oijp^) in the chiral expansion with 
baryon chiral Lagrangian, we can write 

a{n,x,ii) = — + Il{n,x, fi), (72) 
U{n,x,fi) = —y-lJin+^^n (73) 

where x = np/n, / ~ -Ftt and T^kn is the KN sigma term 

^KN^\{m + ms){N\us + ss\N) (74) 

where m = {mu + md)/'2, with the subscripts u, d and s standing, respectively, for up quark, 
down quark and strange quark. In chiral perturbation theory, the third term - known as 
Weinberg- Tomozawa (WT) term - corresponds to 0(p) term and the sigma term to 0{p'^). 
From the point of view of renormalization-group equation, what governs the flow to kaon 
condensation is the WT term with the sigma term being irrelevant. 

Note that the kaon condensation threshold density ra^ for a given ji will be given by 

a{n = rij^, x, fi) = 0. (75) 

We should inject one remark here, returning to more discussions later. In the ab- 
sence of other model-independent tools, chiral perturbation approach has been used to 
address the problem both for finite nuclei and infinite matter, in particular kaon conden- 
sation at high density. Unfortunately lacking lattice information, the reliability of chiral 
perturbation theory in medium, while reasonably successful for free-space processes, has not 
been established. Dense matter is a highly correlated system that is difficult to access by 
perturbative methods. For instance in kaon nuclear systems, while low-order chiral pertur- 
bation calculations, unitarizcd or otherwise, tend to give a shallow kaon-nuclear potential 
whereas phenomenology indicates an attraction several times typical chiral perturbation 
predictions [56]. This deeply attractive potential in nuclear matter can be understood in 
a Walecka-type mean-field approach [57]. Indeed the deeply bound kaon-light nuclear sys- 
tems proposed by Akaishi and Yamazaki [58] require certain nonperturbative mechanisms 
not caught in low-order chiral perturbation expansion. In this connection, there is a hint 
from explorative lattice QCD simulation of kaon condensation [59] on dense kaon systems 
containing up to 12 K~^s which indicate that the properties of the condensate are remark- 
ably well reproduced by leading-order chiral perturbation theory. While direct baryonic 
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background, difficult to implement in lattice calculations, is still needed, what this surpris- 
ing result indicates, somewhat persuasively, is that there can be substantial cancelations in 
(numerous) higher chiral-order terms that cannot be captured by a few perturbatively cal- 
culable terms. This also questions the validity of "higher-order" calculations in discussing 
the depth of kaon-nuclear potentials. 

In the next subsection, we will present one possible non-perturbative effect connected 
with the half-skyrmion phase. 

To proceed, let us simply assume that higher chiral-order and higher-K field terms 
can be ignored, and compute the proton and neutron chemical potentials modified by the 
kaon condensation 



JJ_ ^KN 



P P 



K^, (76) 

(77) 



where 



|u°-^4 = 4{l-2^)SN{nN). (78) 

UN 

In beta-equilibrium with kaon condensation, the baryon number density is still carried by 
the nucleon, n = ns, while strangeness is carried entirely by kaons. If one ignores the 
back-reaction of kaon-nuclear interactions on nuclear interactions, then we can write the 
chemical potential difference using S{p) and K, 

Hn-fip = ^l-2x)SN{n) + ^K\ (79) 

Equation (79) indicates that the symmetry energy will be modified when kaons con- 
dense. To see how it should be modified, we redefine the symmetry energy in the presence 
of kaon condensation as 

E{n,x,K'^) = E{n,l/2,K^) + E,ym{n,x,K^) (80) 
Esym{n,x,K^) = E{n,x,K^)-E{n,l/2,K^). (81) 

Then the symmetry energy density is given by 

€symin,x,K^) = nil - 2xfSN{n) + [//^(n,^) - //^(n, l/2)]i^2, (82) 

where the first term is the symmetry energy of pure neutron matter 

e%^{n,x) = n{l-2xfSN{n) (83) 

and the chemical potential /i(ri, x) is given by the solution of 

^i'-m],+ ^-^^^n + ^n = Q. (84) 
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Figure 5: jjL{n,x) = lok{x.,p) with a parameter set used in [11]. 



The numerical calculations in a model adopted in this work, shown in Fig. 5, can be ap- 
proximated as 

IJL{n,x) ^ a{l — hx)mKi (85) 

where a{n) and h{n) (of order of 10~^ for n/no = 1 — 6) are density-dependent constants. 

On the other hand the kaon chemical potential can be formally derived from the 
Lagrangian we are using as a function of density given by 

M(n,x) = -[-^ (86) 

= a*{l-h*x) (87) 

where uk is the kaon number density and 

* = lf!!^__!Li h* = —i\—-—\ 

The numerical calculations confirm the formal expression and give a* ~ arriK and b* ~ b. 
Then we get 

fi\n,x) - i?{n, 1/2) ~ (1 - 2x)^ml. (89) 

It shows that the x-dependence is not in the form of (1 — 2x)'^ but linear in (1 — 2x), 
which breaks n-p permutation symmetry. The basic reason is that when K~ condenses, 
the isospin is spontaneously broken and the Weinberg-Tomozawa term in Eq. (73) brings 
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in isospin asymmetry. Thus the K condensation gives an additional contribution to the 
symmetry energy in the form 



Esymin,x,K^) = (1 - 2xfSNin) + (1 - 2x)SKin,K), 



(90) 



where 



SK{n,K) 



a?{n)h{n) K'^m\ 
2 ~ 



K 



(91) 



Note that for the parameters given in Fig. 5, Sxin, K) is of order of a few MeV. Hence the 
effect is not significant when compared to Sn. 

Ahhough we do not know exactly how symmetry energy will be modified in kaon- 
condensed matter, one of the major effects of kaon condensation inferred from this simple 
example is that the x— dependence of symmetry energy deviates from (1 — 2x)^. Whether 
it can be measured at heavy ion machine including FAIR, NIC A or KoRIA (or "RAON") 
is an interesting question, if it is plausible that the time scale for nuclear matter at the 
collider is long enough to realize the kaon condensed matter via weak interaction. 

5.2.1 In-medium kaon mass in the skyrmion crystal 

In this subsection, we turn to the skyrmion description and present a nontrivial 
consequence of the half-skyrmion phase on the in- medium mass m'^{= fx) given by 



and on kaon condensation in dense matter. The self-energy 11 is given by (73) to 0{p^) in 
chiral perturbation theory, but as mentioned, it is unlikely that ChPT is reliable when kaons 
can condense. As an illustration of how nonperturbative effects can figure in the process, we 
examine how kaon spectrum is affected by the skyrmion-half-skyrmion transition at density 



(50) as in the Callan-Klebanov bound-state model for hyperons [60, 61]. At high density, 
the skyrmion background is expected to be reliably described by SU (2) skyrmions put on 
crystal as in [62]. We briefly describe what was found in that work and then discuss its 
ramification for the dropping kaon mass in compact-star medium . 

The Lagrangian (50) is extended to three flavors by supplementing with the Wess- 
Zumino term and taking the chiral U field as 




(92) 



ni/2- We consider kaon fields fiuctuating on top of a dense skyrmion background given by 




(93) 




(94) 
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Here uq G SU{2) is the skyrmion background with winding number B. When expanded in 
terms of the kaon field K, the Lagrangian we are concerned with is given by 

^ = ^SV{2) + ■Ck-SU(2) (95) 

where the kaonic part of the Lagrangian (50) in the background of the skyrmion background 
uo(x) and the dilaton xoi ^k-su(2)) is of the form 

'Ck-SU(2) = {D^K)^D^'K - K\ml + V)K + z{K^D^'KB^ + h.c.) + • • • (96) 

where D^K = (9^ + V^)ir with = ^{dfj,uoul + d^uluo) , V is a time-independent potential 
term, is the baryon current and 2; is a known constant in front of the topological Wess- 
Zumino term. Explicitly 

'Ck-su(2) = aid^K^d^K) + iPiK^k - K^K) - jK^K + ■■■ (97) 

where 

a = ^(^^0 + ^^0 + 2), 

P = ^B%uo + ul + 2), 

7 = —mj^{uo + ul + 2) (98) 

where k = Xo/fx the baryon charge density. 

The Lagrangian (97) describes the property of the kaon fluctuating in the background 
of the skyrmion uq and the dilaton xo- We can consider the skyrmion with any baryon 
number A. For A = B = I, hyperons are obtained [60, 61] when collective-quantized by the 
SU{2) rotator coordinate A(t) 

U{r,t) = A{t)Uo{r)A-\t), 

K{r,t) = A{t)U{r,t) (99) 

The mass difference between the proton and the hyperons Am > 175 MeV can be accounted 
for by the binding of up to three kaons. The binding is predominantly given by the Wess- 
Zumino term with the hyperfine splittings obtained from collective-quantization. In terms 
of the Nc counting, the Wess-Zumino term gives the vibrational correction of 0{N^) and 
the hyperfine splitting of 0{N~^) to the baryon mass. In [60], the pure Skyrme model 
was studied with some success. With the introduction of the vector mesons p and u, one 
obtains an even nicer fit to all hyperons in the baryon octet multiplet [61]. It is likely that 
one would obtain an even better description with the infinite tower as argued above but 
this has not been confirmed yet. 

Here we shall discuss applying the above strategy to many-nucleon systems as was 
done in [62] . Details have not been worked out but it is clear that the procedure is straight- 
forward, provided that the uq for multi-baryon system is known. We shall take uq studied 
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in [62]. Ideally it would be best within our framework to embed kaons into A-body sys- 
tems with A — )• oo described by uq simulated with the skyrmions put on crystal using 
the Lagrangian (40). However this has not been done yet. We will therefore consider the 
simpler Lagrangian (50). What we are particularly interested in is the effect of the non- 
perturbative phenomenon associated with the change-over from a skyrmion matter to a 
half-skyrmion matter at ni/2 on the properties of the fluctuating kaons on top of the given 
dense background. 

To have a qualitative idea, one can simplify the calculation by averaging the quantities 
(98) over the density-dependent skyrmion background as well as the dilaton condensate xo- 
Then the s-wave kaon dispersion formula from (97) is 

aojj^ + 2l5ojK + 7 = (100) 

where a, ^ and 7 are the averaged quantities. The solutions are given in Fig. 6 for two 
values of the dilaton mass with typical values for the parameters indicated in the figure 
caption. 




Figure 6: w-^i vs. n/ng (where ng — 0.16 fm^"^ is the normal nuclear matter density) in 
dense skyrmion matter which consists of three phases: (I) {qq) 7^ and /* / 0, (II) {qq) = 
and f^y^O and (III) {qq) = and /* = 0. The parameters are fixed at V^g/tt = rup = 780 
MeV and dilaton mass = 600 MeV ((a) left panel) and = 720 MeV ((b) right pannel). 
It should be noted that while n-^ is very sensitive to m^^, ni/2 is practically independent of 

One can see that the kaon mass drops similarly to the tree-order chiral perturbation 
results (73) up to density ng, then undergoes a propitious drop not present in chiral per- 
turbation theory at due to the appearance of the half-skyrmion phase. Although the 
location of ni/2 is quite insensitive to the dilaton mass as long as it is low ^ 1 GeV, 
the magnitude of the drop however does depend strongly on it. Such a drop - if it is real - 
would represent a genuinely nonperturbative effect brought in by the topology change^^. 

^^This is just about the same amount of a additional attraction required by Akaishi and Yamazaki [58] 
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5.2.2 Hyperons cannot appear before kaon condensation 

Kaons will condense in compact star matter when the electron chemical potential 
fie is equal to the effective kaon energy or mass m*^^ (for s-wave). As discussed above - 
and one can see from Fig. 6, kaon condensation in the skyrmion-crystal description will 
depend sensitively on the precise value of the dilaton mass. Without dilatons, mj^ will 
not fall low enough to condense up to density n ~ 5no, but with the dilaton of ~ 700 
MeV appropriate for nuclear physics [31], < 200 Mev for n < 2.5nQ. This suggests 
that kaon condensation at a reasonable density is likely to take place in this picture. It is 
intriguing that this property in the half-skyrmion phase is quite close to what one finds by 
fluctuating from the vector manifestation fixed point of hidden local symmetry theory as 
discussed in [63].^^ Without the dilaton and the change-over from the skyrmion phase to 
the half-skyrmion phase, this would not occur. This again reminds one of the proximity 
of the half-skyrmion phase to the Georgi vector symmetry connected to the onset of the 
mended symmetry arrived at in the dilaton limit discussed above. 

It should be emphasized that there has been an extensive debate with widely divergent 
conclusions as to what happens to kaon condensation when hyperons start figuring in the 
system. The issue is whether the role of the strangeness hyperons bring in is additional 
to that of condensed kaons and if so, their presence would not delay or even banish kaon 
condensation. The debate on this matter is typically based on mean-field treatments of 
phenomenological Lagrangians or low-order chiral Lagrangians supplemented with "natural" 
higher-order meson field operators. But the problem is that such treatments cannot be 
trusted with confidence. That is because such treatments suffer from one well-known defect: 
One does not know how to write an effective Lagrangian with higher-field operators that give 
in the mean-field a reliable description at higher density than that of normal matter. The 
customary procedure is to mix two unrelated ingredients, one anchored on low-order chiral 
perturbation theory and the other a phenomenologically implemented set of massive vector 
and scalar degrees of freedom. The two ingredients are unconstrained from each other, and 
hence although experiments guide theorists up to rig, there is nothing that guides them 
beyond it, particularly at high density relevant to compact stars. This accounts for the 
widely varying predictions of EoS for n > uq. 

Here we would like to suggest that the Callan-Klebanov bound-state treatment used 
above offers a possible solution to this problem. 

For this purpose, we can continue from Sec. 3.3.7 where the symmetry energy was 
calculated by quantizing the neutron matter in a skyrmion-crystal description. We take uq 

for their compressed nuclei. 

^*The behavior of at n > ni/2 can be compared with the tree-order chiral perturbation with BR 
scaling taken into account - with « 250 MeV, which is consistent with current lattice results (see [11] for 
comments on this matter.). This would imply for electron chemical potential « 130 MeV a critical density 
~ 3.2no, consistent with what was found in [63, 57]. A similar result was obtained by Westerberg [64] 
using the Callan-Klebanov bound-state model embedded in a hypersphere. He found ~ 3.7no. 
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to be a multi-skyrmion configuration for infinite matter put on crystal. We then embed 
kaons into that system and quantize the system using the cohective coordinates (99). In 
this way, binding Nk kaons to the soliton with the winding number ^4 ^ 1 would be 
tantamount to having hyperons in the system with the strangeness Nk- The binding comes 
mostly from the Wess-Zumino term which is the leading-order term involving strangeness, 
going as 0{N'^) in the large Nc limit. To this order, there is no distinction between the 
different hyperons for the given Nk and there are no interactions between them. They will 
however appear at the next order as hyperfine effects accounting for splitting between 
members of equal strangeness. 

In the CK scheme described above, it is easy to see that hyperons cannot appear 
before kaons condense. To see this, it suffices to first recall that kaons condense when the 
in-medium mass of the kaon is equal to the electron chemical potential /Xg 

mk = l^e (101) 

and then note that in the CK model, hyperons appear when kaons get bound to the 
skyrmion, mainly by the Wess-Zumino term. The mass difference between the nucleon 
and the hyperon (denoted Y) in the CK model is 

my-mN^ rnl^ + ©(l/ATj (102) 

where rn^^ is the effective mass of the kaon (< m/^) in the soliton background which is of 
0{N^. ) in the Nc counting. The hyperfine effect is in the 0{1/Nc) terms. In compact stars, 
hyperons will appear when 

my — rriN = IJ-e- (103) 

This is the same condition as kaon condensation (101) since in medium, m^^^ = m'^. Thus 
to 0{N^), hyperons appear at the density at which kaons condense. To go to higher order 
in 1/Nc in medium, one would have to consider loop corrections that are density-dependent 
in addition to the hyperfine effects coming from the collective quantization. This would 
require a lot more work going beyond the present scheme, but it seems likely that the 1/Nc 
corrections will not upset this conclusion. 

5.2.3 Hyperons in standard RMF approaches 

As argued above, in the unified approach based on skyrmion crystal, the role of 
kaon condensation and that of hyperons are intricately related. In this case, it would be 
a double-counting if both were treated independently and combined. In the framework of 
phenomenological approach as well as in relativistic mean field approach, however, such a 
unified treatment of the two related mechanisms does not exist. What one does is to mix 
the two in an arbitrary way, unconstrained by each other. 
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Here we will simply ignore kaons and bring strangeness degrees of freedom only 
through hyperons. We may consider this as a formulation in terms of baryonic degrees 
of freedom only with the kaon degrees of freedom "integrated out." 

For a system of nucleons and hyperons, the hyper charge Y = S + B (or strangeness 
number S) is conserved in strong interactions as are the proton (/a = 1/2) and neutron 
(I3 = —1/2) numbers in SU{2). We can classify the eigenstate by the definite number of 
protons and neutrons and total number of hypercharge (73 = l/2{Np — Nn),Y = Sjl^ATy.): 

\Np,Nn,Y). (104) 

Although the SU{3) symmetry is broken explicitly, the isospin symmetry can be used to de- 
compose the eigenstate, Eq.(104), into the irreducible representations (multiplets) of SU{2) 
as given by 

\Np,Nn,) = J:(^i,Y)C(i,Y)\I,y ■■ Np,Nn) (105) 

where = /(/+ 1). The energy of the eigenstate of the Hamiltonian does not depend on 
Is but on and Y as given by 

E{Np,Nn,Y) = {Np,Nn,Y\H\Np,Nn,Y) = S(j,y)|Cg,y)|%,y) (106) 

where Ef^j^y) is an reduced matrix element of the Hamiltonian, H = Hq + Hi^t, 

E^i^Y) = {I,Y\\H\\I,Y) (107) 

which is independent of Is, and depends on the details of the strong interactions for each 
{I, y)-channel. 

Therefore one can expect the nuclear symmetry energy can be affected by the presence 
of hyperons. The simplest parametrization inferred from the pure nuclear matter is 

E{n, X, Vi) = E{n, x = 1/2, y^) + EgyminN, x, yi), (108) 

where yi are the density of hypcron i and njv is the nuclcon number density njv = Ji-n + ^p- 
n is the total baryon number density n = tin + y, where y = T,iyi = n — tin is the total 
hyperon number density. 

Let us consider first, as a simple example, the non-relativistic and non-interacting 
system. Although the kinetic contribution to the symmetry energy is small compared with 
the interaction term, particularly at high density, one can still gain a valuable insight from 
it into how to parameterize the symmetry energy and how the presence of hyperons can 
affect it. 

Recall that with only two flavors, n and p, the symmetry energy was defined in section 

2.3 as 

E{n,x) = E{n,l/2) + Esym{n,x) (109) 
Esym{n,x) = E{n,x)- E{n,l/2), (110) 
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where E is the energy per baryon. In the non-interacting Fermi-gas description, 

<-(n,x) = Ep[2'/'({l-xf/' + x'/')-l] (111) 



where x = This can be approximated for the full range of x from to 1 in a simple 
form -within 5% accuracy - as 

= Sf^^%n)il-2xf, (114) 

S^-^(n) = (22/3 _ 1)3^0 (^)2/3_ (115) 

5 no 

However when we take the above form with (1 — Ix)"^ as an expansion around x = 1/2, then 
it cannot be naively used near x ~ (pure neutron matter) or 1 for pure proton matter. 
Nonetheless in the literature the symmetry energy of the above form has been employed in 
many of parametrization schemes even for the pure neutron matter. A recent discussion of 
this form for x = — 1 can be found in [16] [17]. We would like to see what comes out for 
the hyperonic matter in this simple free-gas approximation. In fact, we know very little of 
at what density the hyperonic matter can become relevant but let us simply assume that 
we have three types of hadrons, n, p and a hyperon (it can be any of the hyperons for the 
moment). The n-p symmetric state is defined by evenly distributed densities , n„ = n^, in 
the presence of the hyperons with density ny. Then the energy per particle can be written 
as 

E{n,x,y) = E(n,x = 1/2, y) + Esymin,x,y), (116) 
Esym{n,x,y) = E{n,x,y) - E{n,x = 1/2, y) . (117) 

Here y is the fraction of hyperon number density, y = ny/n. In non-relativistic approxima- 
tion, 

nE{n,x,y) = nN^E^Kl - af^ + {1 + af/^] + nyl-^{37rWf/^ (118) 
z zmy 

E{n,x,y) = '^\lE%C^fl\l-afl^ + {l + afl^] + y'l^\^^^^ 

n 2b riQ 5 my no 

= [(l-y)5/3l[(i_„)5/3 + (i+^)5/3]+2V3i!Ly5/3] 3^0(^)2/3 (^20) 

2 my 5 no 

(121) 

with a = 1 — 2x. Since n^v = n — ny = (1 — y)n, the pure nuclear matter would correspond 
to y = 0. The symmetry energy from the kinetic energy part is 

Esym{n,x,y) = E{n,x,y) - E{n,x = 1/2, y) (122) 
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= [(1 - y)V3 1 [(1 _ afl^ + (1 + a)V3 _ 2]]\E%{^fl^ (123) 
z 5 no 

= {l-yf'^Sf'^^{n){l-2xf, (124) 

where Sf'^^^{n) = (2^/^ _ i^3£'0^(^2L)2/3 Thus the effect of hyperons on the kinetic energy 
part of the symmetry energy is to reduce its strength by a factor of (1 — y)^^^. With nuclear 
interactions turned on, the contribution from the interactions will of course change the 
symmetry energy significantly. 

It is convenient to parameterize the density dependence of 5 in a polynomial form as 

SiriN) = ^M—r, (125) 
no 

where is a density- independent constant. In this form, the kinetic energy contribution 
corresponds to k = 2/3, 

S2/3 = sli:' = i2y'~l)lEl (126) 

To estimate the contributions from the nucleon-nucleon interactions, we assume that the 
interactions between nucleons are little affected by the presence of hyperons even beyond 
the hyperon threshold. Then 

Esymin,x,y) = E%^{nN)C^) = E^y^{nN){l-y) (127) 
= {l-2xf~S{n,y). (128) 

This shows that the density in the argument of symmetry energy inferred from pure nuclcon 
matter, E^^^, should be the nucleon number density, not the total baryon number density 
and that the symmetry energy per baryon becomes smaller by a factor of (1 — y) since the 
baryon number is also carried by the hyperons. Now using 

KmM = (1 - 2x)^5(nAr) = (1 - 2x)2s«5«(-)«(l - y)(«+i) (129) 

77-0 

we get 

Sin,y) = S«5«(-)«(l-y)('^+i). (130) 
no 

For the kinetic contribution, we get 

S'^'^'in, y) = {l- yf/^Sf'-^^n) (131) 
as discussed above. What we learn from this simplified calculations is that 

S{n,y)^S{n). (132) 
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5.3 Quark degrees of freedom 

Similarly to the case of strange hadrons in dense hadronic matter, we might wonder 
whether one cannot expect more massive hadronic degree of freedom when density is much 
higher. However, it appears that the next heavy hadrons, i.e., charmed hadrons, are too 
heavy to be excited at the density regime relevant to compact star interior. Instead we 
expect a new form of matter in terms of quarks if the deconfinment phase transition can 
take place before other more massive hadronic degree of freedom appear. The criterion 
for the onset of the phase transition to a quark matter is whether the pressure of hadronic 
matter at a given density is greater than the pressure of quark matter. To address this issue, 
we would need reliable information on quark EoS. The critical density will be determined 
at the phase boundary where the pressures are balanced. In our approach, we will suppose 
this phase change occurs at the end of kaon condensed nuclear matter. 

The density at which quark matter presumably becomes relevant is much higher than 
the normal nuclear matter density, but it is not in the density regime where perturbative 
QCD can be applied. There can be a plethora of domains where hadronic and quarkish 
phases overlap, for which we have neither theoretical control nor experimental data. For 
our purpose, what is needed is a quarkish equivalent to the hadronic symmetry energy. As 
mentioned, the symmetry energy in the hadronic sector is controlled by strongly correlated 
interactions, with the kinetic energy term playing a minor role at high density. Whether 
similar correlations are involved in the quark sector is unknown. For example, there have 
been no serious discussions on the density dependence of the Clebsch-Gordon coefficients, 
C/, for a given number of u, d and s quarks. In the SQM, it is assumed that J = is 
dominated. The quantities related to quark-quark correlations such as A^d or A^j. analogous 
to Ann or Anp has not been sufficiently explored except for the case of the free quark-gas 
approximation. 

We will therefore explore a scenario of the phase change which bypasses the details 
of the EoS of quark matter at the phase boundary with the assumed kaon condensed nuclear 
matter. 

5.3.1 Non-strange quark matter 

Quark matter can appear as a result of confinement-deconfinement phase transition. 
This phase transition is basically a strong interaction process. Suppose nuclear matter with 
no strangeness changes over to a deconfined quark phase. There will be no room for strange 
quarks in this case because it is the strong interaction that is involved. At a later stage, 
it can evolve into a strange quark matter via weak transition emitting neutrons out of the 
system if SQM is more stable energetically for a given quark number density. 

We first consider nuclear matter without kaon condensation (NM for short) and de- 
confined non-strange quark matter (NSQM for short). Suppose there is a phase boundary 
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between NM and NSQM at the density Uc (for NM) and (for QM)-^^, respectively, in 
the interior of a neutron star. The chemical equilibrium reads 

IJ-n = ^IJ-d + fJ'U, IJ-p = fJ'd + '^fJ'U- (133) 
The symmetry energies in each phase are related to the chemical potentials as 

Hn-I^p = 4(1 - 2a;jv)S'jv(njv) (134) 
^^d-^^u = i{l-2xQ)SQ{nQ), (135) 

where xjv = Up/riN, xq = n^/ng and ng = riu + Ud- Prom the chemical equilibrium 

fJ-n-IJ'p = IJ-d-IJ'u, (136) 
using XQ = Hu/uQ = {1 + xn)/S, we get 

(1 - 2xN)SNinc) = ^(1 - 2xN)SQin^) (137) 

from which we obtain a constraint on the symmetry energy at the chiral restoration density 
He (which we assume to coincide with the deconfinement density), 

SN{nc) = Isgin'^) (138) 

for xn 7^ 1/2. Just for an illustration, suppose we take an expression for the symmetry 
energy factor of the non-interacting quark gas: 

5Q = ^(37r^nQ/2)V3[2V3_i]^ (I39) 

and assume that the critical quark number density, ra^/3, is equal to the critical baryon 
number density, ric, then we will get 

SN{nc) = 24.6(— )V3MeV . (140) 
no 

We may obtain a different constraint on the possible form of the nuclear symmetry energy, 
if nuclear matter undergoes a phase transition into a quark phase at pc- However, that 
will yield a highly nontrivial constraint [5]. It also depends on the nontrivial quark-quark 
interactions in the deconfined phase. In a recent paper by G. Pagliara and J. Schaffner- 
Bielich [66], it is argued that the Sq can be three times larger than that given by the free 
quark model if the quark phase is in 2SC, which is of course a highly correlated system, 
unlike free quarks. In Eq. (137), the simplest solution is x;v = 1/2 irrespective of the specific 
form of symmetry energy. We think this possibility can be realized in stellar matter at high 
density. We shall obtain this result from a different consideration. However in heavy ion 
collision, x is fixed to be less than 1/2, under which the phase transition occurs. So we have 
to solve nontrivial phase boundary conditions [5] [6] [66]. 

We are taking into consideration that there may be a density jump between the NM and the QM at the 
boundary, such that ^ Sric- 
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5.3.2 Strange quark matter 

In dense compact-star matter, the system evolves in weak equilibrium, so strange 
quarks figure in dense matter through weak interactions. The nature of interactions for 
the phase change is, however, different if it is from a phase in which kaons induced by the 
weak interaction arc present to start with. Now if kaon condensation has taken place in 
the system, then there is already non-vanishing strangeness in the nuclear matter up to 
the phase boundary. In this case, we can imagine the confinement-deconfinement phase 
transition taking place constrained by the weak equilibrium in the kaon condensed matter, 
leading to strange quark matter. In this subsection, we develop the scenario in which kaon 
condensation joins smoothly to a strange quark matter. 

Consider the case where strange quark matter (SQM for short) meets kaon-condensed 
nuclear matter (KNM for short). As was done in the previous section, we can obtain the 
kaon threshold density, n^, by the profile of the kaon chemical potential. Prom up to 
the critical density for chiral restoration ric, the weak equilibrium condition reads 

li = 4{l-2x)SN{nc) + F{K,^i,nc). (141) 

With the Weinberg-Tomozawa term for the kaon-nucleon interactions introduced in Section 
4.2, F{K,fi,nc) can be expressed as 

F{K,ii,n^) = ixF{K,n^), (142) 

and we have 

/X = 4(1 - 2x)SN{nc) + liF{K, rie). (143) 

As the kaon chemical potential - equivalcntly effective mass - ^ approaches at the critical 
density, the solution x = ^ = 1/2 appears naturally at the phase boundary, provided 
F{K,nc) 7^ 1 which is assumed to be valid for the range of density we are concerned with. 
At the phase boundary, the chemical equilibrium (via confinement-deconfinement) 

reads 

l^n- fJ-p = l^d- l^u, (144) 
I^K- = l^s- l^u- (145) 

Note that the strange quark is required at the boundary, which implies that there should 
be a strange quark matter for n > ric- Since fi{= ixk) = , we have from Eqs. (144) and 
(145) 

liu = l^d and fis = l^u, (146) 

which gives 

IJ'u = li-d = l^s- (147) 
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This is the chemical potential relation for the SQM in the massless limit. In this simple 
picture, the kaon condensed nuclear matter leads naturally to a strange quark matter. In 
this sense, it is natural to consider the kaon condensed state as the doorway to strange quark 
matter as proposed elsewhere. The coincidence of kaon condensation and chiral symmetry 
restoration in the chiral limit was predicted a long time ago with a kaon bound to a skyrmion 
put on hyper sphere [70]. 

6 Discussions and Conclusion 

In this article reviewing our work done in the WCUIII/Hanyang program, we dis- 
cussed how flavor symmetries figure in the nuclear symmetry energy which plays an impor- 
tant role in the structure of neutron-rich nuclei and in the high-density EoS of compact stars 
. In this way of looking at the problem, the prominent role of the tensor forces in nuclear 
interactions is highlighted. When dense baryonic matter is described by multi-skyrmions 
put on crystal, the structure of the tensor forces is controlled by a change in the parameters 
of our generalized hidden local symmetry Lagrangian which is engendered by a change of 
topology from the skyrmion matter to a half-skyrmion matter at a density ni/2 lying slightly 
above the nuclear matter density. We suggest that this density regime could be explored 
by such RIB accelerators as the "KoRIA" (alias "RAON" ) being built in Korea. 

The new scaling of the parameters, called "BLPR scaling" replacing the old BR 
scaling, if correct, will be of crucial importance in the effort to access the regime of density 
where the chiral phase transition or deconfinement of quarks is expected to take place. We 
have found that the new scaling makes the EoS - which is soft below ni/2 - become stiff 
above, precisely the mechanism needed to understand the 1. 97 Mq star without invoking 
other degrees of freedom than nucleonic. On the other hand, the topology change at ni^2 
makes the effective anti-kaon (e.g., K~) mass in medium to undergo a propitious drop 
so that it could speed up kaon condensation, which would then soften the EoS. We have 
also argued that hyperons can appear only when kaons condense, so they need not be 
considered as two different phenomena [72] . Since the stiffening of the nucleonic EoS defers 
kaon condensation to higher density, it is not clear what would be the net consequence 
between these two opposing effects. What we have found however is that even if kaon 
condensation takes place, it can play the role of a "doorway state" to strange quark matter 
and compact star masses of ~ 2Mq can be accommodated. 

We close with some additional remarks. 

One of the observations made with the skyrmion-half-skyrmion transition is the cusp 
in the symmetry energy at ni/2- The same cusp structure is expected when the crystal 
result is transcribed into a many-body theory language in terms of the BLPR and one 
assumes that the symmetry energy is dominated by the tensor forces. Such a cusp was 

^"The relation (147) will of course be spoiled in nature with the s-quark mass. But going smoothly over 
from a kaon-condensed state to strange quark matter should remain intact. 
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not observed in standard phcnomcnological or relativistic mean field treatments. It is 
also possible that such a cusp structure will be smeared out by many-body correlations 
in microscopic many-body or field-theoretic approaches. Such a renormalization-group- 
implemented field theoretic calculation using Viowk with the BLPR suitably incorporated, 
discussed in [50], does indicate the topology change taking place at n^/ is smoothed by 
correlations. 

Another observation that merits to be highlighted is that most of the parametric 
forms we have looked at are found to show A„p, which is negative at normal nuclear 
density, becoming positive at higher density, say, ~ Suq. This implies that the attraction 
between neutron-proton at lower density is changing into repulsion at higher density whereas 
the repulsion between neutron-neutron (A„„) keeps growing. 

An intriguing consequence of the suppression of the p tensor force at density n 1> 77-1/2 
predicted in this paper is that at high density, > 3no, the pion tensor will dominate. This 
could lead to a pion-condensed crystalline spiral structure predicted a long-time ago by 
Pandharipandc and Smith [68]. Given that this structure arises due to the fractioning of 
skyrmions into half-skyrmions, this may be related to the "baryonic popcorn" structure 
discussed in holographic approach to baryons [69] . 

Among the characteristics of the hidden local symmetry approach, a distinctive fea- 
ture that is yet to be confirmed by either theory or experiments is the vector manifestation 
which predicts that at the chiral transition (in the chiral limit), a set of vector mesons 
become massless with a flavor gauge symmetry "emerging" from collective phenomena in 
dense and/or hot medium [23]. Although there are indirect evidences for it in the form 
of "dropping mass," there is no direct or "smoking-gun" evidence either for or against it. 
The highly heralded dilepton production in heavy-ion collisions is not a simple and trustful 
probe for in-medium vector-meson mass - and hence chiral symmetry - contrary to what 
has been naively thought [71]. It remains to find an experiment that will do the job. 
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